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The rotor-router model on a graph describes a discrete-time walk accompanied by 
the deterministic evolution of configurations of rotors randomly placed on vertices of 
the graph. We prove the following property: if at some moment of time, the rotors 
form a closed clockwise contour on the planar graph, then the clockwise rotations 
of rotors generate a walk which enters into the contour at some vertex v, performs 
a number of steps inside the contour so that the contour formed by rotors becomes 
anti-clockwise, and then leaves the contour at the same vertex v. This property 
generalizes the previously proved theorem for the case when the rotor conhguration 
inside the contour is a cycle-rooted spanning tree, and all rotors inside the contonr 
perform a full rotation. We use the proven property for an analysis of the sub- 
diffusive behavior of the rotor-router walk. 
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INTRODUCTION 
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The Eulerian walkers model introduced in [![ as an example of self-organized criticality 
,l3|], was rediscovered later by researchers in different fields and attracted much attention 
due to its simple algorithmic structure 0-ll|. Some properties of the Eulerian walkers, 
stated in [lland in the related paper 12|, were formulated then as rigorous graph-theoretical 
theorems 13|. In the mathematical literature, the model received the name “rotor-router 
walk” proposed by Jim Propp who considered it as a derandomized analogue of the 
ordinary random walk. 

A formal dehnition of the rotor-router walk is given in Section 11. In a less formal way, 
one considers a lattice with arrows attached to the lattice sites. The arrow at every site is 
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directed to one of its neighbors on the lattice. A particle (chip) performs a walk jnmping 
from a site to a neighboring site. Arriving to a given site, the chip changes the direction of 
the arrow at this site to next position, for instance 90 degrees clockwise, and moves to the 
neighbor pointed by new position of the arrow. Thus, given an initial orientation of arrows 
on the whole lattice, the rotor-router walk is fully deterministic. Fig JT] illustrates three steps 
of the rotor walk on the square lattice. 
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FIG. 1: Circles denote the lattice sites, (a) The chip is originally in the filled circle where the 
arrow is directed ”up”. (b) The chip rotates the arrow clockwise and moves right, (c) The next 
clockwise rotation sends the chip down, (d) The last position of the chip is the right lower corner. 


At present, the basic properties of the rotor-router walks on hnite graphs can be con¬ 
sidered as well established. The walk starting from an arbitrary vertex of a hnite Eulerian 
graph G (see the dehnition below) eventually settles into an Eulerian circuit where each 
directed edge of G is visited exactly once. When the walker is in the Eulerian circuit, con- 
hgurations of rotors associated to vertices of G are recurrent. A striking property of the 
recurrent state is a hlling of loops formed by rotors on a planar graph G under clockwise 
rotations of rotors: if a sequence of rotors makes up a closed clockwise oriented contour G 
on G at some moment of the Eulerian circuit, the walker enters into contour G in a vertex 
n G G, performs a full rotation of each rotor internal to G, changes the orientation of rotors 


on G to opposite and leaves the contour in the same vertex v 
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. Il3l. This property of the 

n 

contour is called reversibility. Recently, Chan, Church and Grochow [1^ have noticed that 
the reversibility of all cycles on a graph is a test whether the given graph is planar or not. 
Specihcally, they proved that a connected loopless graph is planar if and only if all cycles 
are reversible under the rotor-router dynamics. 

The behavior of the rotor-router walk on inhnite graphs is quite different. Let G be an 
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infinite square lattice with the rotors associated to each vertex. Each rotor has four possible 
directions ordered clockwise. The initial orientations of rotors are independently distributed 
with equal probabilities for each direction. The walker starting at the origin performs a 
diffusion-like motion. However, the exponent v in the law of the mean-square displacement 
for n steps {r‘^{n)) ~ differs considerably from its random walk value u = 1/2. Extensive 
computer experiments show that z/ = 1/3 with high precision. A qualitative derivation of 
the subdiffusion law z/ = 1/3 in |l| is far from a rigorous proof although the arguments 
given there seem to be credible. These arguments are based on an assumption that the 
walker’s trajectory covers a kernel of the cluster almost densely before reaching the shell of 
the cluster again. Then, the radius of kernel grows with a velocity inversely proportional to 
the area of the kernel, so the average radius R ~ Despite a strong numerical evidence 
for the assumption, there are no theoretical explanations of this unusual behavior of the 
rotor-router walk. 

The aim of the present work is to consider the motion of the rotor-router walk in a random 
environment in more detail. Analyzing the trajectories of a walk inside closed contours, we 
hnd a rigorous rule of visiting a closed contour emerged in a random rotor conhguration. 
Namely, we prove a property which we call the weak reversibility for the rotor-router walks, 
where the condition of the full rotation of each rotor internal to the contour is released. 
Instead, some rotors inside C perform a partial rotation or do not move at all. Nevertheless, 
two properties of reversibility are still retained: the rotor-router walk entering C in a vertex 
V E C leaves the contour in the same vertex v, and the clockwise orientation of rotors on 
the contour C becomes anti-clockwise. 

We consider the rotor-router walk inside the contour and we describe precisely which part 
of the contour interior is hlled by the trajectory of the walk. Then, we discuss a possible 
relation of the loop reversibility to the subdiffusion of the rotor-router walk on the inhnite 
lattice. We observe that the set of vertices, where rotors make up clockwise contours, grows 
according to a dehnite rule, namely, the sequence of such vertices generates a spiral-like 
structure. Moreover, the obtained spirals, being random, in average obey an unexpected 
law which we call the asymptotical Archimedean property. The spiral-like structure of the 
reference points of the closed clockwise contours together with the weak reversibility of 
contours provides the subdiffusion behavior of the rotor-router walk. 
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II. THE MODEL AND A BASIC THEOREM 

Consider a directed graph (digraph) G = (V, E) with a set of vertices V = V(G) and a set 
of directed edges E = E(G). We assume that there are no self-loops or multiple edges in G 
although the dehnition of the model can be extended to this case. If for each edge directed 
from V to w there exists an edge directed from w to v, we call the graph G bidirected. A 
bidirected graph can be obtained by replacing each edge of an undirected graph with a pair 
of directed edges, one in each direction. A spanning subgraph G' of a bidirected graph G is 
a digraph with the set of vertices V(G') = V(G) and a set of edges E{G') C E{G). 

A path of length n from vertex oGCtofoGCisa sequence of distinct vertices 
ui, ^ 2 , • • • Wn+i such that Vi and Uj+i are connected by an edge e* G E, i = l,2,...,n, 
Vi = a, Vn+i = b. The path becomes a cycle if a = 6. A shortest possible cycle has length 
2 and consists of two adjacent vertices Vi, ^ 2 , which are connected by a pair of edges from 
Vi to V 2 and back. We call such cycles dimers by analogy with lattice dimers covering two 
neighboring vertices. A cycle formed by more than two edges is called contour. 

An Eulerian circuit on a hnite digraph is a walk which starts and ends on the same 
vertex and visits each directed edge exactly once. If such a walk exists, the digraph is called 
Eulerian. A digraph is strongly connected if for any two distinct vertices v, w there are 
paths from v to w and from w to v. A strongly connected digraph G = {V,E) is Eulerian if 
and only if for each vertex v G V in-degree and out-degree of v are equal. In particular, the 
one-component bidirected graph is Eulerian. 

The rotor-router model is dehned as follows. Consider an arbitrary connected digraph 
G = (y,E). Denote the number of outgoing edges (out-degree) from the vertex n G E by 
d^. The total number of edges of G is \E\ = Each vertex n G E is associated with 

a rotor, which is directed along one of the outgoing edges from v. The rotor directions at 
the vertex v are specihed by an integer variable which takes values from 0 < q;„ < — 1 

for dv > 1. 

The set p = {Q;t,| n G E, 0 < < dy — 1} dehnes the rotor configuration. Starting with 

an arbitrary rotor conhguration, one drops a chip to a vertex of G chosen at random. At 
each time step the chip arriving at a vertex v, first changes the rotor direction from ay to 
ay + 1 , and then moves one step along the new rotor direction from v to the corresponding 
neighboring vertex. The periodicity of ay is assumed (ay ± d„ = ay). 
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The rotor configuration p can be considered as a spanning subgraph of G (p C G) with 
the set of vertices V{p) = V{G) and the set of directed edges E{p) C E{G) coinciding with 
the rotors. The state of the system at any moment of time is given by the pair (p, v) of the 
rotor configuration p and the position of the chip v ^ V. 

A vertex n G id is called sink if its out-degree = 0. In the absence of sinks, i.e. when 
each vertex has at least one outgoing edge, the motion of the chip does not stop. If iterating 
the rotor-router operation from the state (p, v) eventually leads back to (p, v) we say that 


(p, n) is recurrent; transient otherwise. According to arguments in |l|, the rotor-router walk, 
started from an arbitrary initial state (p, v) on a finite graph, passes transient states and 
enters into a recurrent state continuing the motion in the limiting cycle which is the Eulerian 
circuit of the graph. 


In [l5(] a useful notion of unicycle is introduced (see also [l6j). A connected spanning 
subgraph of a digraph G, in which every vertex has one outgoing edge contains exactly one 
cycle. The state (p, n) is called unicycle if the set of edges E{p) contains a unique directed 
cycle and v lies on this cycle. Then, two basic properties of the rotor-router model on the 
Eulerian graphs can be formulated in terms of unicycles. 

Property A ([l^. Theorem 3.8). Let G be a strongly connected digraph. Then a single¬ 
chip-rotor state (p, v) on G is recurrent if and only if it is a unicycle. 

The rotor states that are not unicycles, are transient. In contrast to recurrent states, 
they appear at the initial stage of evolution up to the moment when the system enters into 
the Eulerian circuit. 

Property B ( 13|, Lemma 4.9). Let G be an Eulerian digraph with m edges. Let (p, n) 
be a unicycle in G. If one iterates the rotor-router operation m times starting from (p, n), 
the chip traverses an Euler tour of G, each rotor makes one full turn, and the state of the 
system returns to (p, n). 

A theorem on reversibility of loops at the recurrent state 3, 13] mentioned in Introduc¬ 
tion reads : 

Theorem 1 ( 12] and [^, Corollary 4.11). Let G be a bidirected planar graph with the 


outgoing edges at each vertex ordered clockwise. Let (p, v) be a unicycle on G with the cycle 
G oriented clockwise. After the rotor-router walk makes some number of steps, each rotor 
internal to G has performed a full rotation, each rotor external to G has not moved, and 
each rotor on G has performed a partial rotation so that G is now oriented anti-clockwise. 
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Our aim in the next section is a proof of the theorem on weak reversibility. 


III. CONTOURS IN RANDOM ROTOR ENVIRONMENT 

The theorem on weak reversibility that we are going to prove, is related to planar digraphs. 
However, we start with a lemma which is valid for general directed planar and non-planar 
graphs. 

Consider a bidirected contour C = (ui, ^ 2 , • • •,'^n) in a digraph G = {V,E), that is a 
contour in which the vertices Vi and Uj+i are connected by two edges ef, e^_^i G E one in each 
direction, I < i < n, where the periodicity of the indices is assumed = u, = ef, 
= e~). Given the rotor-router model dehned on G, we say that the bidirected contour G 
obeys the domino ordering if for each rotor at n*, 1 < z < n there exists a direction a*, such 
that the rotor of. points from Vi to Vi-i and a*. -|- 1 points from n, to Uj+i. The directions 
..., are called negative with respect to G, whereas the directions -1-1,..., -|-1 

are called positive, correspondingly (FigI2]). 



FIG. 2: Bidirected contour obeying domino ordering. 

Lemma 1. Given an arbitrary hnite Eulerian digraph G, let G = (ui,... ,n„) be a bidi¬ 
rected contour obeying domino ordering. Let the rotor-router walk starts at the vertex Vn 
from an initial rotor conhguration with positive directions of all rotors at Uj (z = 1, 2,..., n). 
Then, after some number of steps, the walk produces a conhguration with negative direc- 
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tions ..., The moments ti {i = 1, 2,..., n), when the directions a*, are reached, 
are ordered as follows: 0 < < tn-i < • • • < ^2 < < |-E|. 

Proof. Since the digraph G is hnite, each vertex n G C is visited by the walk indehnitely 
many times. Let T be the hrst moment of time when the rotor at returns to the initial 
direction + 1. None of the edges of G is visited twice by the walk before T. Indeed, 
assume the hrst repeated visit happens for an edge e E E directed from a vertex v ^ Vn- 
Then, vertex v is visited more than times, i.e. more than once from one of its neighbors 
V 7 ^ u, that contradicts to the assumption that the hrst edge, visited twice, is the edge 
directed from vertex v. Each rotor being in direction a* + 1 at some moment of time t will 
change its direction to a* if and only if the vertex v is visited by the walk exactly — 1 
times. The hrst vertex v of the contour G being visited dy — 1 times, is the vertex Vn, since 
other vertices u* of G are not visited neither from Uj+i nor from Uj_i before, and there are no 
repeated visits from any vertices. The moment when the direction of arrow at becomes 
is tn < T. The next vertex is Vn-i, as it has been already visited from Direction 
is reached at moment f„_i,where tn < tn-i < T because the arrow at Vn cannot return 
to the initial direction a*^ + 1 earlier than the arrow at vi is in direction a*^. The process 
continues for vertices Vn- 2 , Vn- 3 ,... ,Vi. Non of rotors at u*, z = 1,..., u has direction a* + 1 
before ti < T < \E\, since the vertex Vi is not yet visited from Uj_i □. 

If (p, v) is a unicycle on the planar graph G with the contour G, then the internal subgraph 
Gint C G formed by rotors inside C is a spanning forest, i.e. the graph whose vertices coincide 
with internal vertices of G and edges form trees rooted at G. Now, consider the situation 
when the rotors at the internal vertices of G have arbitrary orientations. Since each internal 
vertex contains a rotor, the internal subgraph G^t remains a spanning subgraph. The edges 
of Gint can be grouped into disjoint cycles and trees rooted at these cycles. The trees which 
are not rooted at the internal cycles, have roots at the external contour G. 

All contours can be either oriented clockwise or anti-clockwise. By an analogy with the 
unicycle, we introduce the multicycle as a graph containing exactly k cycles together with 
k chips at vertices oq, oi,..., a^-i belonging to the cycles. For multicycles, we will use the 
notation (p, uq, ..., 

First, we consider a situation where the external contour is clockwise and all internal 
contours are anti-clockwise (Figj3]). We prove a theorem on weak reversibility generalizing 
the Theorem 1. 


a 
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FIG. 3: Multicycle: external clockwise contour containing the vertex a and three anti-clockwise 
contours containing vertices b, c, d. 

Theorem 2. Let G be a connected bidirected planar graph and (p, oq, ai,..., Ok-i) be a 
mnlticycle with the external contour Cq oriented clockwise together with k — 1 internal cycles 
Gi,..., Gfc_i oriented anti-clockwise. The rotor-router operation is sequentially applied to 
the chip at Oq G Gq until the moment when the chip returns to Oq, and the rotor at Oq is 
made oriented anticlockwise. Then, the same is applied to chips at oi,..., a^-i until the 
moments when chips starting from G Gj return to Oj and the rotors at Oj are made oriented 
clockwise. Then, all rotors on Gq are becoming oriented anticlockwise, whereas all rotors 
on Gi,..., Ck-i become oriented clockwise, also all vertices internal to Go and external to 
Gi,..., Gfc_i perform a full rotation. 

Proof. Given a multicycle (p, Oq, Oi,..., Ok-i), we construct an auxiliary graph G* per¬ 
forming the following operations (FigHj): 

1. reverse all the clockwise rotors on the external contour Gq to anticlockwise, also all 
anticlockwise rotors on internal cycles Ci,, Ck-i to clockwise; 

2. remove all the vertices and edges of G external to Go, also remove the vertices and 
edges internal to any of Gj, i = 1,..., k — 1; 

3. add an additional contour Gq consisting of k vertices Oq, ..., cik-v 

4. connect vertices Oi and a\ by edges for alH = 0,..., fc — 1; 

5. change the directions of rotors at ao,ai,... ,ak-i so that new directions are toward 
Oq, a'l ,..., the order of directions of rotors at Oq, %,..., is shown in FigJUfor the 
vertex o'q. 

6. put a single chip at Oq. 

The obtained single-chip-and-rotor state on the auxiliary graph G* is the unicycle U*. 


9 


After these preparations, we can use the property B for the unicycle U*. According to this 
property, the chip starting at Og returns to Og, whereas the unicycle U* returns to its initial 
state. The initial state of rotors in U* coincides with the hnal state of rotors at vertices of 
Co,..., Ck-i and also at vertices internal to Co and external to Ci,..., Ck-i claimed in the 
theorem. 

Thus, we have only to prove that, during their evolution, these rotors reach the initial 
state indicated in the theorem. The additional contour Cg provides the walker successive 
visits of the vertices Oq, Oi,..., ak-i from Og, ..., ak-v straightforward to see that 
the first visit of the rotor walk to Co is followed by |Co| steps which make the orientation of 
rotors clockwise, whereas the hrst visit to every internal contour C* is followed by \Ci\ steps 
which make the orientation on Cj anticlockwise. Since the sequences of |Co| steps for the 
contour Co and |Cj| steps for contours Cj both are deterministic, we can ignore all of them, 
replacing each sequence by reversed orientations of contours just before the first visit to Co 
and Cj,i = l,...,/c — 1. Thus, the conditions of the theorem are fulhlled. □ 

Corollary. Let ui, ^ 2 , • • •, '^n be vertices of Co ordered clockwise, and Vn = no be the first 
vertex of Co visited by the rotor-router walk. Then the moments of time ti,... ,tn, when 
the rotors at Vi,V 2 , ■ ■ ■ ,Vn become anticlockwise with respect to Cq, are ordered as follows: 

tn ^ tn—1 ^2 ti- 

Proof. After operations 1 and 2, listed in the proof of Theorem 2, the bidirected contour 
Co obeys domino ordering. Then, the statement of corollary follows from Lemma 1. □ 

The condition of anticlockwise orientation of the internal cycles Ci,..., Ck-i in Theorem 
2 can be released. Indeed, let Cj, 1 < j < n — 1 be a cycle oriented clockwise. If there 
are no internal cycles inside Cj, then the hrst hit of a chip to the contour Cj converts, by 
Theorem 1, the orientation of Cj into anticlockwise after visiting its interior. Then, we can 
simply skip the interval of the chip evolution between the moments of entry and exit from 
Cj considering the orientation of Cj as anticlockwise. If the interior of cycle Cj contains 
anticlockwise cycles only, we can skip the mentioned interval of the evolution, by Theorem 2. 
If there are clockwise cycles inside Cj, then we can apply the skipping procedure recursively. 
Since each generation of cycles is enclosed in the previous one, we eventually reach a smallest 
cycle which does not contain any cycle inside. 
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FIG. 4: Unicycle constructed from the external clockwise contour and three anti-clockwise contours. 
Numbers 1,2,3 at Oq correspond to the order of bonds adjacent to this vertex. 

IV. WEAK REVERSIBILITY AND THE SUBDIFFUSION LAW 

As mentioned in the introduction, the crucial concepts used for explanation of the subd¬ 
iffusion behavior of the rotor-router walk on the inhnite lattice are the kernel-shell structure 
of the cluster of visited vertices and the assumption that the trajectory of the walk covers 
the kernel almost densely in intervals between periodic returns to the shell. This property 
was approved by extensive simulations [l| , but yet a detailed explanation has not been pre¬ 
sented. In this section, we use Theorem 2 on weak reversibility to show how the rotor-router 
dynamics leads to the subdiffusion law. 

Let G = (U, E) be an inhnite square lattice and po be the initial rotor conhguration 
with {a^\v G U} taken uniformly from the set {0,1,2,3}. The rotor-router walk starts the 
motion from the origin and performs T steps forming a cluster of visited vertices and edges. 
We hx all moments of time when the rotors form clockwise contours, and numerate these 
contours as they appear. Assume that after gi^gpg^ rotor-chip conhguration creates 
a clockwise contour Ck- The T^b"-bth step is directed to the vertex Vk G Ck- According to 
Theorem 2, the chip returns to Vk after visiting the interior of Ck, and the clockwise contour 
becomes anticlockwise. We denote by the moment of exit from Ck and put a label Sk 

at Vk- Continuing, we obtain a sequence of contours Ci,C 2 , ■ ■ ■ and labels si, S 2 ,.... We are 
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interested in relative disposition of contours and labels appeared during time T. 

To distinguish between rotors which did not move until the given moment of time and 
those involved into the motion, we call the moved rotors activated. A current position of 
all activated rotors, except the last one, forms a single rooted tree T. The location of the 
root is a current position of the chip. If the position coincides with the label Sk, the root 
Sk of the tree T is connected with the origin by a radial branch 0 —)■ s^- For t > one 

possibility is that the diffusing chip reaches a vertex connected with the branch 0 —)■ by 
a sequence of activated or non-activated rotors so that the next appeared clockwise contour 
is adjacent to Ck- Otherwise, the chip creates one or several isolated clockwise contours 
having no common edge with the branch 0 —)■ one by one. When, eventually, the chip 

reaches a vertex connected with the branch 0 —?• s^, the next clockwise contour adjacent to 
Ck appears, and all isolated contours (which have already changed their orientation from 
clockwise to anti-clockwise by this time) become connected with Ck- 

Consider a label Sk situated near the boundary of the cluster of activated rotors. A 
preferable position for s^+i is on the right side of the branch 0 —)■ sa, to provide clockwise 
orientation of the contour Ck+i if it has common edges with 0 —)■ s^. Then, the preferable 
direction of successive positions of labels Sk, Sk+i, Sk+ 2 , • • • is clockwise with respect to the 
origin of the cluster. Since the size of cluster grows with time, the positions Sk, Sk+i, Sk+ 2 , ■ ■ ■ 
form a spiral-like structure (Fig. |5]). 

A spiral is called Archimedean if, in polar coordinates {r,6), it can be described by the 
equation 

r = a + b9, (4.1) 

where a and b are real numbers. 

If a spiral is random, we will say that it is asymptotically Archimedean in average if 

(r) 

—)■ 6 for 0 —)■ oo, (4.2) 

where the average is taken over the uniformly distributed states of the spiral. Above defi¬ 
nitions of the sequence of contours Ci,... ,Cn and labels Si, ..., allow us to formulate a 
key conjecture: 

Conjecture 1. Given a random configuration of rotors on the infinite square lattice, 
the spiral of labels si,S2,S3,... generated by the rotor-router walk is asymptotically 
Archimedean in average taken over the uniformly distributed initial states of rotors. 
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FIG. 5: The spiral-like sequence of labels in a cluster obtained after T = 10^ steps of the rotor- 
router walk. Rainbow colors from violet to red correspond to increasing time. 

A numerical verification of the conjecture is not a trivial task. The simulation of 10"^ 
walks of length T = 10® shows that the ratio r{sn)/0{sn) as a function of n for a single 
spiral fluctuates around a constant value up to very large n. However, for the ratio of the 
root-mean-square values, we estimated the convergence to the constant as: 

~ 1.85 + 0(n-2). (4.3) 

^{OisnY) 

Now, we can represent the rotor-router walk as a sequence of transitions between labels 
Si —)■ S 2 —^ • • • —^ Sn and the time intervals [T 2 *”\ ... which the chip 

spends for visiting the interiors A(Gi), A{C 2 ),... of contours Gi, G 2 ,.... For each A: > 1, 
there exists k' such that the walk between and s^' is a loop of the spiral. Consider a 
domain constituted of contours associated with labels s^, s^+i,..., s)., 

k' 

(4,4) 

j=k 


Vi 


The area of is of an order where R is the loop radius. By Theorem 2, the rotor- 
router walk visits interiors of all contours in A^ one by one and, therefore, the total time 
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AT needed for closing the loop is of an order B?. An advance Ai? of the loop radius for the 
time interval AT is of the order of the spiral step 27r6. Thus, the velocity of growth of R is 
proportional to the inverse area of the loop, 


AT 1 


(4.5) 


from which we obtain the asymptotical law R{T) 2^1/3 fQj, large T. 

The arguments leading to (14.51) are qualitative, as well as those in [l|. However, now 
we are equipped with three new concepts, namely: the weakly reversible contours, the 
labels marking the clockwise contours and the spiral-like structure of labels. They make the 
statements of [l| more constructive and show clear meaning of the shell-kernel structure of 
the growing clusters used for the explanation of the subdiffusion law. 

We can conclude that the statistics of labels appears to be a basic feature of the rotor- 
router walk, which determines its long range behavior. Beside the spiral-like structure, the 
labels can be characterized by the average time interval between them. Specihcally, the 
average number of steps = T,[+] — between labels and Sn+i, be. the interval 

between exiting from the clockwise contour and entering to Cn+i, tends to a constant 
value for large n: 

{5tn) ~ 6.81 -I- 0{n~^). (4.6) 


Since the spiral step tends to the constant, the spacial average density of the labels p{r) 
does not depend on distance r from the origin in the large interval of r, 1 <C r -C Our 
Monte-Carlo simulations conhrm this property. We found that near the origin, for which 
the spiral structure is not pronounced yet, the average density has a peak with p(0) ~ 0.37, 
which decays rapidly to the plateau value pst — 0.13. For r > the density vanishes 

sharply in accordance with the subdiffusion law T(T) 2 ^i/3_ 
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